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Abstract 

We investigate the geometry of injection regions and its relationship to optimization of power flows in tree 
networks. The injection region is the set of all vectors of bus power injections that satisfy the network and operation 
constraints. The geometrical object of interest is the set of Pareto-optimal points of the injection region. If the voltage 
magnitudes are fixed, the injection region of a tree network can be written as a linear transformation of the product 
of two-bus injection regions, one for each line in the network. Using this decomposition, we show that under the 
practical condition that the angle difference across each line is not too large, the set of Pareto-optimal points of the 
injection region remains unchanged by taking the convex hull. Moreover, the resulting convexified optimal power 
flow problem can be efficiently solved via semi-definite programming or second order cone relaxations. These results 
improve upon earlier works by removing the assumptions on active power lower bounds. It is also shown that our 
practical angle assumption guarantees two other properties: (i) the uniqueness of the solution of the power flow 
problem, and (ii) the non-negativity of the locational marginal prices. Partial results are presented for the case when 
the voltage magnitudes are not fixed but can lie within certain bounds. 

I. Introduction 

AC optimal power flow (OPF) is a basic problem in power engineering. The problem is to efficiently allocate 
power in the electrical network, under various operation constraints on voltages, flows, thermal dissipation and 
bus powers. For general networks, the OPF problem is known to be non-convex and is challenging U), 12. Some 
of the earlier analysis has focused on understanding the existence and the behavior of load flow around local 
solutions Q, 10. Recently, different convex relaxation techniques have been applied to the OPF problem in an 
attempt to find global solutions 0, J6). It was recently observed in Q that many practical instances of the OPF 
problem can be convexified via a rank relaxation. This observation spurs the question: when can an OPF problem 
be convexifed and solved efficiently? This question was partially answered in several recent independent works 
ll8l- lfT0l : the convexification of OPF is possible if the network has a tree topology and some conditions on the 
bus power constraints hold. The goal of this paper is to provide a unified understanding of these results through a 
deeper investigation of the underlying geometry of the optimization problem. Through this understanding, we are 
also able to strengthen these earlier results. 

There are three reasons why it is worthwhile to focus on tree networks. First, although OPF is traditionally solved 
for transmission networks, there is an increasing interest in optimizing power flows in distribution networks due to 
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the emergence of demand response and distributed generation lUTl. Unlike transmission networks, most distribution 
networks have a tree topology. Second, as will become apparent, assuming a tree topology is a natural simplification 
of the general OPF problem ( |[T2l made a similar observation for the power flow problem), and results for this 
simplified problem will shed light on the general problem. Third, as was shown in (9), if one is allowed to put 
phase shifters in the network, then the OPF problem for general network topologies can essentially be reduced to 
one for tree networks. 

Following (8), our approach to the problem is based on an investigation of the convexity properties of the power 
injection region. The injection region is the set of all vectors of feasible real power injections Pj's (both generations 
and withdraws) at the various buses that satisfy the given network and operation constraints. We are particularly 
interested in the Pareto-front of the injection region; these are the points on the boundary of the region for which 
one cannot decrease any component without increasing another component. The significance of the Pareto-front 
is that the optimal solution of OPF problems with increasing objective functions defined on the injection region 
must lie there. The first question we are after is: although the injection region is non-convex, when does its Pareto- 
front remain unchanged upon taking the convex hull of the injection region? This property would ensure that any 
OPF problem over the injection region is convexifiable: solving it over the larger convex hull of the injection 
region would yield an optimal solution to the original non-convex problem (in the sequel, we will abbreviate this 
property by simply saying that the Pareto-front is convex). While convexifiability is a desirable intrinsic property 
of any optimization problem, a second question of interest is: can the resulting convexified OPF problem be solved 
efficiently? 

To answer these questions, the first step is to view the injection region as a linear transformation of the higher 
dimensional power flow region: this is the set of all vectors of feasible real power flows P^'s, one along each 
direction of each line. We first focus on the case when the voltage magnitudes are fixed at all buses. In the space 
of power flows, the network and operation constraints in general networks can be grouped into three types: 

1) local constraints on the two flows and P^i along each line (i, k): these include angle, line flow and thermal 
constraints. Figure Q] gives an example of the (nonlinear) feasible set of {Pik,Pki) due to flow constraints. 
Note that all these local constraints are effectively angle constraints. Also this feasible set can be interpreted 
as the power flow region of a two-bus network with impedance given by that of the line (i, k). 

2) global constraints on the flows due to bus power constraints. These are linear constraints. 

3) global Kirchhoff constraints on the flows due to cycles. These are non-linear constraints. 



Fig. 1. The feasible set for the two flows along a line when there are power flow constraints. It is a subset of an ellipse which is the feasible 
set when there are no constraints other than the fixed voltage magnitudes at the two buses. In this example, the feasible set is part of the 
Pareto-front of the ellipse. 
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The third type of constraints is most complex; they are global and non-linear. By focusing on tree networks, we 
are left only with constraints of type 1 and type 2. In this case, it is easy to see that the power flows along different 
lines are decoupled, save for the global (but linear) bus power constraints. The overall power flow region is thus 
simply the product of the two-bus power flow regions, one for each line, intersecting the bus power constraints. 

By exploiting this geometric structure of the power flow region, we answer the first question posed above: if 
the two-bus power flow region associated with each line itself has a convex Pareto-front, then the overall injection 
region has a convex Pareto-front. Thus, a local convexity property guarantees a global convexity property. It is 
shown that the local property holds whenever the angle difference along every line is constrained to be not too 
large (say less than 45°). As we will argue, this holds in most cases in practice due to line flow, thermal or explicit 
angle constraints because of stability reasons. 

Concerning the second computational question, we observe that in our geometric picture, a semi-definite program- 
ming relaxation corresponds to taking the convex hull of the purely voltage-constrained injection region followed 
by intersection with the local and global constraints. This convex relaxation results in general in a set larger than 
the convex hull of the injection region. However, it turns out that our analysis of the first question in fact implies 
that the Pareto-front remains unchanged even with this more relaxed convexification. This provides a resolution to 
the computational question. 

The present work improves upon the earlier papers on tree networks ISl- lTTOl in two ways. First, the arguments 
used in those papers are algebraic and some used non-trivial matrix fitting results, while the present paper uses 
entirely elementary geometric arguments. This geometric approach provides much more insight on the roles of the 
various types of constraints and also explains how the assumption of tree topology simplifies the problem. Second, 
the convexity results in all of the earlier papers require some restriction on the bus power lower bounds (no lower 
bounds are allowed in |9j, iflOl . and any two buses that are connected cannot simultaneously have lower bounds in 
[8)). The results in this paper require no such conditions. Instead, they are replaced by constraints on the differences 
between voltage angles at adjacent buses, which we verify to be satisfied in practice. This latter condition imposes 
a local constraint on the power flows, and is discovered through the geometric decomposition of the power flow 
region. 

We also show that the angle assumption gives rise to two other important properties: (i) the solution of a power 
flow problem becomes unique, and (ii) the locational marginal prices (LMPs) can never become negative. The LMP 
is a common pricing signal used in practice for charging customers and paying generators located at various buses. 
It is known that in congested transmission networks LMP at some buses could be negative ff3l . In this undesirable 
situation, some consumers will be paid to withdraw power. We show that this situation does not happen for a tree 
network with realistic angle assumptions. 

The paper is organized as follows. In Section HU we state the physical model used in the paper. Section [III] 
focuses on the case when the voltages magnitudes at all buses are fixed. We first start with a two-bus network with 
angle, thermal and flow constraints. Then, we consider a general tree network with only local constraints. Finally, 
we add the global bus power constraints to arrive at our full results. We also study the implications of our result 
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on the uniqueness of power flow solutions and the non-negativity of locational marginal prices. Section [TVl extends 
some of the results to networks with variable voltage magnitudes. Section [V] concludes the paper. Similar results 
can be derived if the network has reactive power constraints in addition to active power constraints. Due to space 
limitations, we do not include them here. 

Summary of notations we use throughout the paper: 

• vectors and matrices: We use the notations x and X to denote vectors and matrices, respectively. Given 
two real vectors x and y of the same dimension, the notation x < y denotes a component-wise inequality. 
We denote Hermitian transpose of a matrix with (-) H and conjugation with conj(-). The vector x y is the 
component-wise product of the two vectors x and y, and diag(X) returns the vector containing the diagonal 
elements of the matrix X. The notation j is reserved for \f—l in this work. 

• sets: We use scripted capital letters A,B, . . . to represents sets, which are assumed to be subsets of R n unless 
otherwise stated. Given a set A, conv(^4) denotes the convex hull of A. A point x E A is Pareto-optimal if 
there does not exist another point y E A such that y < x with strict inequality in at least one coordinate. Let 
0(A) denote the set of all Pareto-optimal points of A, which is sometimes called the Pareto front of A. Note 
that if a strictly increasing function is minimized over A, its optimal solution must belong to O(A). 

II. Model 

Consider an AC electrical power network with n buses. With no loss of generality, we assume that the network 
is a connected graph. Following the convention in power engineering, complex scalars representing voltage, current 
and power are denoted by capital letters. We write i ~ k if bus i is connected to k, and i k if they are not 
connected. We often regard the network as a graph with the vertex set V = {!,... ,n} and the edge set £. For 
example, the notation (i, k) € £ implies that there exists a line connecting bus i and bus k. Let denote the 
complex impedance of the line (i, k) and yn- = = go- — jbik represent its admittance, where gik,bik > 0. 
Define the admittance matrix Y as 

Yji^i Va + Vii if i = k 



tti~h, (!) 



-Uik 

if i <■*< k 

where ya is the shunt admittance to ground at bus i. Note that this matrix is symmetric. 

Let v = (Vi, V2, ■ ■ ■ , V n ) G C" be the vector of complex bus voltages and i = (ii, 1-2, ■ ■ . , I n ) € C" be the 
vector of complex currents, where Ii is the total current flowing out of bus i to the rest of the network. By Ohm's 
law and Kirchoff's Current Law, i = Yv. The complex power injected to bus i is equal to Si = Pi + jQi = Vilf , 
where and Qi denote the net active and reactive powers at this bus, respectively. Let p = (P 1; P 2 , . . . , P n ) be the 
vector of real powers, which can be written as p = Rc(v conj(i)) = Re(v (Y H v H )) — Rc(diag(vv ff Y ff )). 



March 2, 2013 



DRAFT 



5 



III. Fixed Voltage Pareto Optimal Points 
A. Two-Bus Network With Angle, Thermal and Flow Constraints 

Consider the two-bus network in Figure [2] with the line admittance g — jb. Let the complex voltages at buses 

\i g-jb V 2 



Fig. 2. A two-bus network. 



1 and 2 be expressed as V\ = \Vi \ exp(j'#i) and Vi — | V2 1 cxpQ'^)- Throughout this subsection, assume that the 
magnitudes |Vi| and are fixed, while 9\ and 82 are variable. The power injections at the two buses are given 
by 



Pi = \Vi\ 2 g+\Vi\\V 2 \bsm(9) - |T4||y 2 | 5 cos(^) 
Pi = \V 2 \ 2 9- \Vi\\V a \ban(9) - IVMgcostf), 



(2a) 
(2b) 



where 8 = 81—62- Since the network has only two buses, P\ = P\ 2 and P2 = P21, where Pj& is the power 
flowing out of bus i to bus k. Since the voltage magnitudes are fixed, the power flows between the buses can both 
be described in terms of the single parameter 8. Notice that a circle centered at the origin and of radius 1 can be 
parameterized as (cos(0), sin(0)). Therefore, d2} represents an affine transformation of a circle, which leads to an 
ellipse. This ellipse contains all points (Pi,/^) satisfying the inequality 



b -9 


-1 r 


-b -g _ 





\Vi\\V 2 \, 



Pi-\Vi\ 2 g 
P2-\V 2 \ 2 g 

where j| • H2 denotes the 2-norm operator. As can be seen from the above relation, the ellipse is centered at 
(|yi| 2 <7, |V2| 2 <7), where its major axis is at an angle of —45° to the x-axis with length IV1V2I& and its minor 
principle axis has length | Vi V2 1 3 ■ If the line is lossy, the injection region is a hollow ellipse as shown in Figure 
|3(a)| If the line is lossless, the ellipse is degenerate and collapses into a line through the origin as shown in Figure 
|3(b)| In practice most lines in distribution networks are lossy with b/g ratio typically between 3 to 5 (instead of 





(a) lossy 



(b) lossless 



Fig. 3. The region defined by I0: (a) shows the region corresponding to |Vi| = | V2 1 = 1 (per unit), b = 5 and g = 1; (b) shows the region 
for a lossless line. 

> 10 in transmission networks) lfl4l . ifTSI . Thus, the interesting and practical case is when the region is a hollow 
ellipse. Note that the convex hull of this region is the filled ellipse. 
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Now, we investigate the effect of thermal, line flow and angle constraints. Since the network has fixed voltage 
magnitudes, the thermal loss and line flow constraints can be recast as angle constraints of the form 9_< 8 < 8 for 
some limits 9 € [— ir, 0] and 8 <G [0, ir]. More precisely, the loss of the line, denoted by Li 2 , can be calculated as 

L 12 = \Vi-V 2 \ 2 g = P12 + P21 

(3) 

= \V 1 \ 2 g-2\V 1 \\V 2 \gsYn{8) + \V 2 \ 2 g. 
It follows from the above equality that a loss constraint L\ 2 < L\ 2 (for a given Li 2 ) can be translated into an 
angle constraint. Likewise, the line flow inequalities P\ 2 < P\ 2 and P 2 i < P 2 i are also angle constraints. As a 
result, we restrict our attention only to angle constraints in the rest of this part. 

We define the injection region V to be the set of all points {(Pi,P 2 )} given by © by varying 8 € [8_,8]. The 
bold curve in Figure [TJ represents the injection region after a certain angle constraint. 

The key property of the non-convex feasible set V for a two-bus network is that it has the same Pareto-Front 
as its convex hull (see Figure [T). To understand the usefulness of this property in solving an optimization problem 
over this region, consider the following pair of optimization problems for a strictly increasing function /: 

minimize /(Pi, P 2 ) (4a) 

subject to (Pi, P 2 ) e V, (4b) 

and 

minimize /(Pi, P 2 ) (5a) 

subject to (Pi, P 2 ) 6 conv(P). (5b) 

Since / is strictly increasing in both of its arguments, the optimal solution to (O must be on the Pareto boundary 
of the feasible set; therefore both optimization problems share the same solution (P* ,P 2 *) <E V. This implies that 
instead of solving the non-covex problem (|4), one can equivalently solve the optimization © that is always convex 
for a convex function /. Hence, even though V is not convex, optimization over V and conv(P) is equivalent for 
a broad range of optimization problems due to the following lemma. 

Lemma 1. Let V £ M 2 be the two-bus injection region defined in (f2|i by varying 8 over [8_, 9]. The relation 
0(V) = C(conv(P)) holds. 

B. General Network With Local Constraints 

In this subsection, we extend Lemma Q] to an arbitrary tree network with local constraints while section IIII-DI 
states and proves the general result with both local and global constraints. 

First, we express the injection region of a general tree as a linear transformation of the power flow region. Given 
a general network described by its admittance matrix Y, consider a connected pair of buses i and k. Let P^ denote 
the power flowing from bus i to bus k through the line (i, k) and P/., denote the power flowing from bus k to bus 



March 2, 2013 



DRAFT 



7 



i. Similar to the two-bus case studied earlier, one can write: 

Pik = \Vi\ 2 g ik + \Vi\\V k \bik sin 6 ik - \Vi\\V k \gik cos 6 ik 

Pki = \V k \ 2 g l k - \Vi\\Vk\bik sin 9ik - |ViI|Vfc[ff<fc cos9 ik , 

where Oik = Oi — 9k- The tuple (Pik, Pki) is referred to as the flow on the line (i,k). As in the two-bus case, all 
the thermal and line flow constraints can be cast as a constraint on the angle 9ik. Note that the angle constraint on 
9ik only affects the flow on the line (i, k); therefore it is called a local constraint. 

There are 2\£\ numbers describing the flows in the network. Let T denote the feasible set of the flows in IR 2 I £ ', 
where the bus voltage magnitudes are fixed across the network and each flow satisfies its local constraints. Recall 
that the net injection at bus i is related to the line flows through the relation Pi = J2 k -k~i ^ ik - This motivates the 
introduction of an n x 2\£\ matrix A defined below with rows indexed by the buses and the columns indexed by 
the lines: 

!\ if { — k 
(7) 
otherwise. 

The matrix A can be seen as a generalization of the edge-to-node adjacency matrix of the graph. The injection 
vector p G R™ and the flow vector f G T are related by p = Af. We express the set of line flows as {Pik,P k i} 
and say that p is achieved by the set of flows. This implies that the feasible injection region V is given by 

V = AT. (8) 

Since the above mapping is linear, it is straightforward to show that conv('P) = Aconv(J r ). 

We now demonstrate that T has a very simple structure: it is simply a product of the two-bus flow regions, one 
for each line in the network: 

t= n (9) 

(i,k)eS 

where the two-dimensional set Tik is the two-bus flow region of the line (i, k). In other words, the flows along 
different lines are decoupled. To substantiate this fact, it suffices to show that the flow on an arbitrary line of the 
network can be adjusted without affecting the flows on other lines. To this end, consider the line (i, k) and a set 
of voltages with the angles 9\, . . . , 9 n . The power flow along the line (i, k) is a function of = 9i — 9 k , Assume 
that we want to achieve a new flow on the line associated with some angle 9 ik . In light of the tree structure of the 
network, it is possible to find a new set of angles 9i , . . . , 9 n such that 9, L — 9 k = 9 ik and that the angle difference 
is preserved for every line in £\(i, k). 

Due to this product structure of F, it is possible to generalize Lemma [TJ 

Lemma 2. Given a tree network with fixed voltage magnitudes and local angle constraints, consider the injection 
set V defined in ©. The relation 0(V) = 0(conv(V)) holds. 

Proof: First, we show that 0(conv('P)) C 0{V). Given p G 0(conv('P)), let {(P ik ,Pki)} € conv(J") be the 
set of flows that achieves p. Consider a line (i, k) G £. Since J 7 is a product space and p G 0(conv('P)), we 
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Pi 



(a) (b) 
Fig. 4. Three possible cases for the bus power constrained injection region. 



(O 



have (Pik,Pki) & O (cony (IFik)). Moreover, it follows from Lemma [T]that O(conv(j-ifc)) = O(Fik)- Therefore, 
(Pik, Pki) S Fik for every line (i, k). This gives p e V and consequently p € 0(V). 

Next, we show that 0(V) C C(conv(P)). Given p e O('P), assume that p ^ (9(conv('P)). Then, there exists 
a point p' S 0(conv('P)) such that p' < p with strict inequality in at least one coordinate. By the first part of the 
proof, we have p' £ V, which contradicts p G 0(V). ■ 

C. Two-Bus Network with Bus Constraints 

So far, we have studied tree networks with local angle constraints and without global bus power constraints. We 
want to investigate the effect of bus power constraints. We first consider the two-bus network shown in Figure |2] 
and incorporate the angle constraints together with the bus active power constraints of the form P^ < Pi < Pi 
for i = 1, 2. Let Ve = {(Pi, P2) : |Vi| = Vi, |V 2 | = V 2 , < < 9} be the angle-constrained injection region, 
and Vp = {(Pi,p2) '■ P_i < Pi < Pi, i = 1,2} be the bus power constrained region, where V\ and V2 are the 
given nominal values of the voltage magnitudes. The overall injection region is given by the intersection of the two 
regions through the equation 

v = v g riV P . (io) 

There are several possibilities for the shape of V, as visualized in Figures |4(a)| |4(b)| and |4(c)| In Figure |4(a)] both 
buses have power upper bounds. In Figure |4(b)| P\ has upper bound, while P2 has both upper and lower bounds. 
In Figure [4(c)l both buses have lower bounds. It can be observed that 0(V) = 0(comr(V)) for Figures |4(a)|4(b)| 
but this desirable property does not hold for Figure [4(c)| 

Figure [4(c)| means that in the presence of active power lower bounds, the relationship 0(V) = 0(conv(V)) does 
not always hold. This is the reason for the various assumptions made about bus power lower bounds in ll8l- lfT0l . 
Note that 612 in Figure [4(c)1 is allowed to vary from — it to n. However, the angles are often constrained in practice 
by thermal and/or stability conditions. For example, the thermal constraints usually limit the angle difference on a 
line to be less than 10°. Figure [7] in the appendix shows a typical distribution network together with its thermal 
constraints, from which it can be observed that each angle difference is restricted to be less than 7°. Flow constraints 
also limit the angle differences in a similar fashion. 

Assume that the angle constraints are such that Ve = 0(comr(Vg)), implying that every point in Ve is Pareto- 
optimal. Now, there are two possible scenarios for the injection region V as shown in Figure [5] In Figure |5(a)| 
some of the points of the region Ve remain in V and they form the Pareto-front of both V and conv('P). In Figure 
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(a) Feasible (b) Infeasible 

Fig. 5. Either V = 0(P) = 0(conv('P)) or the injection region is empty. 

|5(b)[ V = so then conv('P) = as well. We observe in both cases that 0(V) = 0(conv('P)). Therefore, we 
have 0(V) = 0(conv(V)) if Ve = O(conv(V )). 

In terms of the line parameters bi 2 and 512, the condition Ve = 0(conv(Ve)) can be written as: 

- tan - ^ — ) < 9 12 < 6 12 < tan _1 (— ). (11) 
312 512 

Observe that tSu-aT 1 {b V2 / g V2 ) is equal to 45.0°, 63.4° and 78.6° for ^ equal to 1, 2 and 5, respectively. These 
numbers suggest that the above condition is very practical. Note that since the inductance of a practical AC 
transmission line is larger than its resistance, the above requirement is met if \0 12 \, \9±2\ < 45.0°. It is noteworthy 
that an assumption 

-tan-^f^) <6 l2 <9 12 < tan" 1 ^) (12) 
012 012 

is made in Chapter 15 of Ifl6l , under which a practical optimization can be convexified (after approximating the 
power balance equations). However, our condition (flU is far less restrictive than (Q~2). To understand the reason, 
note that the value is around 0.1 for a typical transmission line at the transmission level of the network |[T6l . 
Now, our condition allows for an angle difference as high as 80° while the condition reported in lfl6l confines the 
angle to 6°. 



D. General Tree Networks 

In this section, we study general tree networks with local angle constraints and global bus power constraints. For 
every bus i € V, let Vi denote the fixed voltage magnitude \Vi\. Given an edge (i, k) £ £, assume that the angle 
difference 6^ belongs to the interval [6 ik , Oik)], where 6 ik £ [— tt,0] and 9ki £ [0,7r]. Define the angle-constrained 
flow region for the line (i, k) as 

F 9lk = {(P ik ,P ki ) : e_ lk < e lk < e ik , m = v u \v k \ = V k } 

The angle-constrained injection region can be expressed as Ve = AFe, where J~e = Ylu k)es -^Sik- Following the 
insight from the last subsection, we make the following practical assumption 

-tan - ^ — ) <6 lk <6 lk <tan™ 1 ( — ), V(i,fc)e£. (13) 
gik gik 

This ensures that all points in the flow region of every line are Pareto optimal. As will be shown later, this assumption 
leads to the invertibility of the mapping from the injection region Ve to the flow region Tq, or equivalently the 
uniqueness of the solution of every power flow problem. 
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Assume that the power injection p must be within the interval [Pj,Pj] for every i £ V. To account for these 
constraints, define the hyper-rectangle Vp = {p : p < p < p}, where p = (P_ 1 , . . . ,P n ) and p = (Pi, . . . , P n ). 
The injection region V is then equal to Ve DVp. In what follows, we present the main result of this section. 

Theorem 3. Suppose that V is a non-empty set. Under the assumption \13\ , the following statements hold: 

1) For every injection vector pg?, there exists a unique flow vector f 6 J such that Af = p. 

2) V = 0{V). 

3) 0(V) = 0(conv(P)). 

In order to prove this theorem, the next lemma is needed. 
Lemma 4. Under the assumptions of Theorem\3\ the relation 0(V) = 0{comi{Ve) n Vp) holds. 

The proof of this lemma is provided in the appendix. Using this lemma, we prove Theorem [3] in the sequel. 

Proof of Part 1: Given peP, consider an arbitrary leaf vertex k. Assume that i is the parent of bus k. Since k 
is a leaf, we have P&j = P&, and subsequently P^ can be uniquely determined using the relation Fjk = 0(Fjk). 
One can continue this procedure for every leaf vertex and then go up the tree to determine the flow along each line 
in every direction. 

Proof of Part 2: Since V is a subset of Ve, it is enough to show that Ve = 0(Ve)- To prove this, the first 
observation is that Te = 0(Fq) = 0(conv(J r e)). Given a point p £ Ve, let f £ Te be the unique flow vector 
such that Af = p. There exist strictly positive numbers {cik, (i, k) £ £} such that f is the optimal solution to the 
following optimization problem 

f = axgmin ^ c *kPik- (14) 
teJr ° (t.k)ee 

Since minimizing a strictly increasing function gives rise to a Pareto point, it is enough to show that there exists 
a set of positive constants c\, Ca, • • • , c n such that the optimal solution of the above optimization does not change 
if its objective function (fl4] i is replaced by X)(ifc)e£ c *^* fc = 2~Zi=i °iPi- Since Te is a product space, we can 
multiply any pair (cik,Cki) by a positive constant, and f still remains an optimal solution. Assume that the tree is 
rooted at 1. Let i be a leaf of the tree and consider the path from 1 to i. Without loss of generality, assume that 
the nodes on the path are labeled as 1, 2, . . . , i. By setting c\ as C12, one can define C2 , . . . , Cj according to the 
following recursion 

Cfe,fe-1 
Cfc-l,fc 

where k ranges from 2 to i. After defining c\, . . . , q, we remove all lines of the path l-i from the network. This 
creates i disconnected subtrees of the network rooted at 1, For each of the subtrees with more than 1 node, 

one can repeat the above cost assignment procedure until c%, . . . , c n have all been constructed. This completes the 
proof. 

Proof of Part 3: For notational simplicity, denote conv(Pg) n Vp as iS. To prove this part, we use the relation 

V C conv(P) C S (15) 
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and the result of Lemma |4] i.e., 

0{V)=0(S) (16) 

The first goal is to show the relation 0(V) Q O(conv(p)) by contradiction. Consider a vector p G 0(P) such 
that p ^ 0(conv('P)). There exists a vector p' G 0(conv('P)) such that p' < p with strict inequality in at least 
one coordinate. Hence, it follows from (fT5t that p is not a Pareto point of S, while it is a Pareto point of P. This 
contradicts ( TToT i. To prove the converse statement 0(conv('P)) C 0(P), consider a point p G 0(conv(P)). In light 
of <Q3), P belongs to S. If p G C(5), then p G 0("P) due to (O. If p <^ 0(5), then there must exist a point 
p' G C(5) = ©(T 3 ) such that p' < p with strict inequality in at least one coordinate. This implies p' G P and 
consequently p' G com^T 5 ), which contradicts p G 0(conv('P)). ■ 



E. Numerical Algorithms for Convexification 

The goal of this part is to understand how the results of the preceding subsection can be used to numerically 
solve an optimization with the feasible set P. The relation 0(P) = 0(coiw(P)) derived before states that the 
minimization of an increasing function over either the non-convex set P or the convexified counterpart conv('P) 
leads to the same solution. However, employing a numerical algorithm to minimize a function directly over conv('P) 
is difficult due to the lack of efficient algebraic representations of conv('P). 

To address this issue, one can decompose P as PgflPp and then use the fact that conv('Pe) and conv("Pp) both 
have simple algebraic representations. Consider the OPF problem 



subject to P_, L < Pi < Pi 
Pi = Y,Pik, 



i G V 
i G V 



(17a) 

(17b) 
(17c) 



{Pik,Pki) G Fe ih , 



(i,k) G £; 



(17d) 



where fi(-) is both monotonically increasing and convex for every i G V, P L and Pi are the bus active power lower 
and upper bounds. Algebraically, the set Fe ik can be represented as all vectors (Pik, Pm) satisfying the relations 



bik —gik 

—bik —Qik 



Pik - \Vi\ 2 g ik 



PM-\V k \ 2 g. 



ik 



m\v k \ 



(18) 



and 



Pki <P ki + ^ ki ^ (Pik -P ik ), 
rik — £_ik 



(P_ik-,P.ki) an d (Pik, Pki) denote the flows associated with line (i, k) G £ for 6ik equal to 6_ ik and 6^, respectively. 
The set Pg is a linear mapping from the product flow region Fg, and hence also has an algebraic representation. 
The convex hull of Fg ik (similarly Tg, Pg) is obtained by changing the equal to sign in ( fT8l to a less than or 
equal to sign. If we replace (I17db with (Pik, Pki) G conv(J r e ifc ), a convex problem is obtained. This convexified 
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(a) conv(-p) (b) corwVe n V P 

Fig. 6. The sets conv('P) and conv'Pe H Vp for a two bus network. The sets are different, but they share the same Pareto Front. 

OPF can be solved as an second order cone problem (SOCP) efficiently in polynomial time Q, (9), ifTTl . Since 
any SOCP problem can be written as a semi-definite programming (SDP) problem, the convexified problem can 
also be interpreted as an SDP |8], l[T8l . 

The feasible region of the convexified problem is not conv('P) = conv("Pg n Vp) but instead comr(Ve) H Vp. 
In general, the convex hull operation and the intersection operation do not commute. However, by Lemma [4] the 
Pareto fronts of conv('P) and co\w{Ve) H Vp are identical (see Fig. |6j. Therefore if the original problem in ( fl2b 
is feasible, optimizing over conv('P) or conv(Ve) n Vp yield the same solution. Moreover, 

« If the solution of the convexified OPF is not a feasible point of OPF, then the original OPF problem is infeasible. 

• If the solution of the convexified OPF is a feasible point of OPF, then it is a globally optimal solution of the 
original OPF problem as well. 

F. Nonnegative Locational Marginal Prices 

The locational marginal price at a given bus (the price per one unit of power) is determined by first adding 
one unit of load to that particular bus and then measuring the increase in the optimal generation cost of the OPF 
problem. The objective of this part is to prove that LMPs are always nonnegative for a tree network under the 
realistic angle assumption cTT~3T > . 

We formalize the definition of LMPs for a nonconvex OPF problem in the sequel. Our definition will be consistent 



, let /.* de- 



with the existing ones for a linearized OPF (such as DC OPF) Q~3). Given a vector e 
note the optimal objective value of the OPF problem dTvT > after perturbing its load vector prj = 
as prj + e. To have a meaningful definition of LMPs, we must assume that /* is differentiable at the point e = 0. 
A set {Ai, A„} is called the set of LMPs for buses 1, 2, .., n if J2iev ^ i6i ^ s ^ e first-order approximation of 
/* — /o . As can be seen from this definition, the set of LMPs is unqiue and well-defined. 

Theorem 5. Under the angle assumption ( 1751 ), the following statements hold for every i £ V: 

1) The LMP Xi is equal to the Lagrange multiplier for the power balance equation P(j i — Pp >i = i Pik m 
the convexified OPF problem. 
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2) The LMP A 2 ; is nonnegative. 

Proof of Part 1: Since the angle constraint (l3[ does not depend on pn, it follows from the argument made in 
Section IIII-EI that /* is equal to the optimal value of the convexified OPF problem after perturbing pu as p^> + e. 
Now, Aj being the Lagrange multiplier for the power balance equation Pc i — Po t = 2~2k~i P* k * s an immediate 
consequence of the well-known sensitivity analysis in convex optimization |fl9l . 

Proof of Part 2: Denote the Lagrange multipliers for the inequalities P_ G . < P Gi and P Gi < P Gi in the convexified 
problem as A € and Xi, respectively, for every i e V. We use the superscript "*" to denote the parameters of the 
(convexified) OPF problem at optimality. It is straightforward to verify that A^ = /'(Pq.) — Ai + A,. To prove the 
theorem by contradiction, assume that some LMPs are strictly negative. In line with the proof of Lemma [4] it can 
be shown that (see d30ll): 

(P*k , Pki ) = arg min A, P lk + \ k P kl 

(P ifc ,P M )econv(.F e .. ) 

(19) 

= arg min XiP lk + X k Pki 

for every (i, k) £ £. Define T as a connected, induced subtree of the network with the maximum number of vertices 
such that A; < for every i £ T. A node k £ V\T is called a neighbor of T if (i, k) £ £ for some i £ T. Define 
T e as the subgraph induced by T and its neighbors. For every line (i, k) £ £, let P ik denote the maximum possible 
flow from node i to node k on this line, i.e., Pi k = max{P ik \(P ik , P ki ) £ J^eik}- F° r every (i, k) £ T e , it can be 
deduced from ( fT9] l and the geometric properties of conv(J r e ifc ) and Fei k (as studied earlier) that 

P* k = P ik if ieT,k#T 

(20) 

p;,. i',: or /^ = p fc4 ifz,fcer 

We orient the edges of T e to obtain a direct graph T e using the following procedure: for every line (i, k) £ T e , the 
orientation of this edge in T e is from vertex i to vertex k if P* k = Pi k . Since T e is a directed tree, it must have a 
node i whose in-degree is zero. Due to (|20K node i must belong to T. Now, one can write: 

P* k = P ik , Vk~i (21) 

On the other hand, since i £ T, we have Xi < and therefore P G . = P G (note that Xi = f'(P G .) — X { + Aj). 
This implies that 

p Gi =J2 T * 

As a result of this equality, if the load Pjj t is perturbed as Prj i + e for a very small number e < 0, the OPF 
problem becomes infinite (because when generator i operates at its minimum output with e = 0, the flows on all 
lines connected to bus i will hit their maximum values). This contradicts the differentiability of /* at e = 0. Hence, 
Ai must be nonnegative. ■ 
The convexification method proposed in J7), J9), ifTol for the OPF problem relies on a load over-satisfaction 
assumption. In the language of our work, this is tantamount to the non-negativity of the LMPs. As a by-product of 
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Theorem [5] we have shown that this load over-satisfaction assumption made in previous papers is satisfied as long 
as a practical angle assumption holds. 

IV. Variable Voltage Pareto Optimal Points 

So far, we have assumed that all complex voltages in the network have fixed magnitudes. In this section, the 
results derived earlier will be extended to the case with variable voltage magnitudes under the assumption P € = — oo 
for every i 6 V. The goal is to study the injection region after imposing the constraints 

Pi<Pi, ieV (22a) 

8ik € l&kJik], (i,k)e£ (22b) 

where \0 ik \, \9i k \ < 90°. Note that the results to be developed next are valid even with explicit line flow constraints. 

Given a bus i £ V, let V_i and V denote the given lower and upper bounds on \Vi\. In vector notation, define 
v = (V_i, . . . , V_ n ) and v = (Vi, . . . , V n ). Given a vector v £ M + , define Vg(v) as the angle-constrained injection 
region in the case when the voltage magnitudes are fixed according to v, i.e. \Vi\ = Vi for i = 1, n. Let V and 
Ve denote the regions for the case with variable voltage magnitudes. One can write V = Ve H Vp, where 

Ve= (J V e (y) (23) 

v<v<v 

The problem of interest is to compute the convex hull of Ve- However, the challenge is that the union operator does 
not commute with the convex hull operator in general (because the union of two convex sets may not be convex). 
In what follows, this issue will be addressed by exploiting the flow decomposition technique introduced in J9). Let 
14.2 denote the convex set of 2 x 2 positive semidefinite Hermitian matrices and Ti. n denote the set of all n x n 
Hermitian matrices. Given a matrix W £ H n together with an edge (i, k) £ £, define: 

• Wik'. (i, k) entry of W. 

• Wife: The 2x2 submatrix of W corresponding to the entries (i,i),(i,k),(k,i),(k,k). The matrix is 
called an edge submatrix of W. 

Define also 

U lk {v) = jw £ % : Wj fc £ Hf, W u = V?, W kk = V k 2 , 
tMiik) x Im(Wi fe ) < Re(VK lfc ), 
Re(W: ifc ) < tan(^fe) x lm{W lk ) 

It can be shown that for every matrix W £ Hifc(v) with the property Rank(W;fc) = 1, there exists an angle 
Oik G \0 ik ,6i k ] such that 

ik 



wue ik 

ViVUOu V k 2 
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Thus, 

J-^(v) = { Re(diag(W tt Y#)) : W e U ik {w), 

Rank(W ifc ) = 1} 

Vik ~Uik j— i 

where = (see J8J, J9J). The flow region .Fe ifc (v) can be naturally convexified by dropping 

—Vik Dik 

its rank constraint. However, the convexified set may not be identical to conv(J r g ifc (v)). We use the notation 
conv(Je it (v)) for the convexified flow region, which is defined as 

conv(^(v)) = { Rc(diag(W 4fc Yi)) : W e H ik (ir)} (24) 
The following sets can also be defined in a natural way: 

coBv(Je(v))= conv(J8, t (v)), 

(i,fc)££ 

conv(7 :, e(v)) = Aconv(J r e(v)) 
Note that conv(-) and conv(-) were the same if the angle constraint (122b) did not exist. 

Lemma 6. Given a vector v, the following relations hold: 

0(Ve(y)) = 0(conv(Pe(v))), (25a) 
0{V{v)) = 0{mm{v e {w)) n V P ), (25b) 
I^J 5mv(?8(v)) = Convex set. (25c) 

v<v<v 

Proof: The proofs provided in Section [HI] for the case with fixed voltage magnitudes can be easily adapted 
to prove (125 al l and ( I25bb . Therefore, we only prove (125c) here. It follows from (l24l that there exists a linear 
transformation from %ikiy) to comr(Fe ik (v)), which is independent of v. We denote this transformation as 
conv(J-e iJt (v)) = lik(rlik{v)). Define "H(v) as f)^ fe s ££ 'Hife(v) and the function /(•) as the natural extension 
of ijfe. Hence, conv(J r e(v)) = l(H(v)). One can write: 

U coEv(P e (v)) = (J Aconv(J- e (v)) 

v<v<v v<v<v 



= A \J W(v)) 

V<V< V 

= Ai( |J W(v) 



(26) 

V<V< V 



v<v< V 

On the other hand, U v< - < -H(v) is a convex set because it consists of all Hermitian matrices W whose entries 
satisfy certain linear and convex constraints. Due to the convexity of this set as well as the linearity of A and I, it 
can be concluded from (|26T i that U v <v<v conv('Pfl(v)) is convex. ■ 
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As pointed out before Lemma [6] conv(-) and conv(-) are equivalent if the angle constraint ( 122b) is ignored. In 
this case, it follows from ( I25cb and the relation 

PjC (J conv("P (v)) C conv(V g ) 

V<V< V 

that conv(Pg) = U v <v<v conv^^v)). In other words, as long as there is no angle constraint, the convex hull 
operator commutes with the union operator when it is applied to (123 i . Motivated by this observation, define conv("Pe) 
as the convex set U v <v<v C0lrv (Ve (v)). We present the main theorem of this section below. 

Theorem 7. For a tree network, 0(V) = 0{co\\v{V)) = 0(comj(V 9 ) n Vp). 

Proof: Since 

V C conv(P) C conv('Pe) H V P , (27) 

it suffices to prove that 0(V) = 0(coTw(Ve) H Vp) (see part 3 of Theorem|3). First, we show that 0(conv(Ve) n 
■Pp) C O(-P). Consider a vector p in 0(conv(:P e )n:Pp). By the definition of convtPe), p e C(conv(P e (v))n7'p) 
for some v. Hence, by Lemma [6] p £ ©(T^v)) and consequently p G V. Now, it follows from d2Tb and p G 
C(convCPe) n Pp) that p e 0{V). The relation 0("P) C O(conv(7' e ) n Pp) can be proved in line with the proof 
of Lemma [4] ■ 

A. Convexification via SOCP and SDP Relaxations 

Theorem E] presents two relations 0(V) = 0(conv(P)) and 0(T) = 0(conv('P £ )) n Vp). Although the first 
relation reveals a convexity property of V, it cannot be used directly to convexify a hard optimization over V. Instead, 
one can deploy the second relation for this purpose. By generalizing the argument made in Subsection IIII-EI we 
will spell out in this sequel how to convexity an OPF problem with variable voltage magnitudes. Consider the 
following OPF problem: 

min 53/i(P G< ) 

iev 

subject to P Gz < P Gt , i £ V 

{Pik,Pu) e ^(v), (i,k)e£ 

V < V < V 

Note that visa variable of this optimization accounting for the voltage magnitudes to be optimized. To convexify this 
optimization, it suffuse to replace the constraint (Pjjt, Pm) £ Fe ik (v) with (Pik, Pu) G conv(J r e ifc (v)). Theorem|7] 
guarantees that OPF and convexified OPF will have the same solution. As shown in Section UlI-EI for a special case, 
the convexified OPF is an SOCP problem. The statement that OPF and its SOCP relaxations own the same global 
solution has been proven in J5] using a purely algebraic (rather than geometric) technique. This result is closely 
related to the prior work E], fl8], iTTO) . which shows that OPF and its SDP relaxation have the same solution. This 
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SDP relaxation can be obtained by dropping a single rank constraint Rank(W) = 1 as opposed to a series of rank 
constraints on the edge submatrices of W. Note that 

• As shown in (9), if the feasible sets of the SDP and SOCP relaxations are projected onto the space for bus 

injections, they will lead to the same reduced feasible set (this property is true only for tree networks). 
« Despite of the fact that these two relaxations convexify the injection region in the same way, it is much easier 
to solve the SOCP relaxation, from a computational point of view. This is due to the total number of variables 
involved in the optimization. 

V. Conclusion 

This paper is concerned with understanding the geometric properties of the injection region of a tree-shaped power 
network. Since this region is characterized by nonlinear equations, a fundamental resource allocation problem, named 
optimal power flow (OPF), becomes nonconvex and hard to solve. The objective of this paper is to show that this 
highly nonconvex region preserves important properties of a convex set and therefore optimizations over this region 
can be cast as convex programs. To this end, we have focused on the Pareto front of the injection region, i.e., the 
set of those points in the injection region that are eligible to be a solution to a typical OPF problem. First, we 
have studied the case when the voltage magnitude of every bus is fixed at its nominal value. Although the injection 
region is still non-convex, we have shown that the Pareto fronts of this set and its convex hull are identical under 
various network constraints as long as a practical angle condition is satisfied. This implies that the injection region 
can be replaced by its convex hull in the OPF problem without changing the global solution. An implication of 
this result is that to convexify the OPF problem, its nonlinear constraints can be replaced by simple linear and 
norm constraints and still a global solution of the original problem will be attained. The injection region of a power 
network with variable voltage magnitudes is also studied. 

Appendix 

A. Thermal Constraints of Distribution Networks 

Every transmission line is associated with a current limit, which restricts the maximum amount of current that 
can flow through the line. Once this limit and the length of the line are known, we can convert it first into a \I\ 2 r 
thermal loss constraint and then into an angle constraint on the line. In what follows, we compute some numbers 
for the 13-bus test feeder system given in lfl4l . This system operates at 2.4 KV line to neutral. The angle limits 
are provided in Figure [7] A pair [a, j3) is assigned to each line in this figure, where a shows the angle between 
the two related buses under typical operating conditions (as given in the data) and j3 shows the limit from thermal 
constraints. 

B. Proof of Lemma [4] 

To prove Lemma [4] we first show that O(S) C 0(V) (recall that S = conv(Ve) fl Vp). Consider a point 
x € 0(S), and denote its corresponding line flow from bus i to bus k with Xik for every (i, k) € £ . Due to the 
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Fig. 7. This figure illustrates the angle constraints in a distribution network. 

relation V = 0(V) derived in Part 2 of Theorem [3] it is enough to prove that x£?. Since x is a Pareto point of 
the convex set S, it is the solution of the following optimization 

n 

x= argmin ^""^ qPj 

pGconvfPe) i=1 

subject to P { < Pi < P l i= 1, 2, n. 
for some positive vector (ci, . . . , c„). To simplify the proof, assume that all entries of this vector are strictly positive 
(the idea to be presented next can be adapted to tackle the case with some zero entries). By the duality theory, 
there exist nonnegative Lagrange multipliers A l5 A n and Ai, X n such that 

n 

x = argmin + Aj - AJP,; - XiPi + A^ 

peconv(Vg) i=l 

or equivalently 

x= argmin I Cj ^ P lk \ (28) 

P econv(p e ) i=1 \ feeV: kr ^ J 

where Cj = c, + A; — A^ for every i e V. By complementary slackness, whenever Cj is less than or equal to zero, 
the multiplier A^ must be strictly positive. Therefore 

Xi = P_i whenever Cj < (29) 

On the other hand, since conv^s) = Aconv(J r 6/) and Fg = Y[(i k)e£ ^ik* ^ results from (l28l that 

(x lkl x ki )= argmin c t P tk + c k P ki (30) 

for every (i, k) € f . In order to prove x <G V, it suffices to show that (xi k ,x k i) € J 7 ^.- Notice that if either Ci > 
or c k > 0, then it can be easily inferred from (l30~t that (xi k ,x k i) £ J~g ik . The challenging part of the proof is 
to show the validity of this relation in the case when Ci,c k < 0. Consider an arbitrary vector y (not necessarily 
distinct from x) belonging to V. Since (y ik ,y k i) € Te ik , it is enough to prove that (x ik ,x ki ) = (y ik , y ki ) whenever 
Ci,c k < 0. This will be shown below. 

Consider an edge (i,k) £ £ such that Ci,c k < 0. There exists at least one connected, induced subtree of the 
network including the edge (i, k) with the property that c r < for every vertex r of this subtree. Among all such 



March 2, 2013 



DRAFT 



19 



subtrees, let Q denote the one with the maximum number of vertices. We define two types of nodes in Q. A node 
r G Q is called a boundary node of Q if either it is connected to some node Z G V\Q or it is a leaf of the tree. We 
also say that a node r G V\Q is a neighbor of V if it is connected to some node in V. By d29] i, if r is a node of 
Q, then y r > P_ r = x r . Without loss of generality, assume that the tree is rooted at a boundary node of Q, namely 
node 1. 

Consider an edge (r, Z) of the subtree Q such that node Z is a leaf of Q and node r is its parent. First, we want 
to prove that yi r > x\ r . To this end, consider two possibilities. If Z is a leaf of the original tree, then the inequality 
(|29] i yields y\ r = y\ > P_ t = x\ — x% r . As the second case, assume that I is not a leaf of the original tree. Let m 
denote a neighbor of Q connected to Z. By analyzing the flow region for the line (Z, m) as depicted in Figure [8(a)] 
it follows from ( f30b and the inequalities c; < 0, c m > that (xi m , x m i) is at the lower right corner of Tq x . Thus, 




(a) (b) 

Fig. 8. Figure (a) shows the flow region for the line (I, m), where (x; m , x m ;) lies at its lower right comer due to c; < and c m > 0. Figure 
(b) shows the flow region for the line (r, I) to illustrate that x r i > y r i (due to T^ rl = 0{comr{Tg rl )) and yi r > xi r ). 

%lm ^ Vmi because of (yi m , y m i) € Fe lm - Let A/j denote the set of all nodes connected to Z that are neighbors of Q. 
One can derive the inequality xi m > yi m for every m G A//. Combining this set of inequalities with x\ = P ; < y\ 
or equivalently 

Xi = ^ x hn + Xlr < ^ Vlm + Vlr = VU 

yields that yi r > xi r . As illustrated in Figure [8(b)] this implies that x r i > Vrl- This line of argument can be pursued 
until node 1 of the tree is reached. In particular, since node 1 is assumed to be a boundary node of Q, it can be 
shown by induction that xu > yu for every node Z such that (1,1) G £. On the other hand, 

^2 xu=x 1 =P 1 <yi= ^2 yu 
lev-. ;~x lev-. i~i 

Therefore, the equality xu = yu must hold for every Z ~ 1. By propagating this equality down the subtree Q, we 
obtain that Xik = yik and xu = yu- This completes the proof of the relation 0(S) C OiJ*). 

In order to complete the proof of the lemma, it remains to show that 0(V) C 0(S). To this end, assume by 
contradiction that there is a point p G ©(T 5 ) such that p ^ 0(S). In light of V C S, there exists p' G 0(S) such 
that p' < p with strict inequality in at least one coordinate. However, since 0(S) C ©(T 5 ), p' belongs to V . This 
contradicts the assumption p G 0(V). 
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